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The recent experimental observation of a metal-insulator transition in two dimensions prompts 
a re-examination of the theory of disordered interacting systems. We argue that the existing theory 
permits the existence of a metallic phase and propose a number of experiments such as magneto- 
conductance and tunnelling in the presence of a parallel field, which should provide diagnostic tests 
as to whether a given experimental system is in fact in this regime. We also comment on a generic 
flow diagram which predicts a maximum metallic resistivity. 
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The discovery by Kravchenko et al. Jl],|| of a metal- 
insulator transition (MIT) in a 2 dimensional system (Si- 
MOSFET) and its confirmation by other workers using 
different device designs || and materials have gen- 
erated much excitement because the conventional wis- 
dom has been that all states are localized in two dimen- 
sions. Up to now the discussion of this phenomenon has 
been based on the scaling theory of localization of non- 
interacting particles, [||f7j even though the possibility of 
unusual superconductivity Q or spin orbit scattering |^] 
has also been raised. On the other hand, within the scal- 
ing theory which includes the combined effect of inter- 
action and disorder |HJ[nj a 2d disordered system may 
remain metallic even in the limit of zero temperature p2| . 
In 2d the expansion parameter is the dimensionless resis- 
tance per square Rn defined as g = ^Ra- For weak ^ 
disorder (g <C 1) the scaling is towards a metallic state 
(dRa/dT > 0) p1JT2|. Furthermore, the theory predicts 
that a magnetic field, via the Zeeman splitting, will drive 
the system towards an insulating state fll| , |l3| . This is in 
agreement with experiment [Q. It is therefore useful to 
revisit this theory in light of the recent experimental de- 
velopment. One reason why the theory has not received 
general acceptance is that the scaling equations have the 
peculiar feature that the scaling variables diverge at some 
finite value of the length scale and the theory becomes 
uncontrolled. While this is certainly true in the vicinity 
of the MIT where g ~ 1 , in this paper we reconsider the 
problem of 2d metallic behavior and argue that for weak 
disorder the theory remains under control over a large 
temperature range, provided the renormalization of the 
energy scale (relative to the length scale) is taken into 
account. In fact this renormalization allows the possi- 
bility of a metallic state with finite resistance in 2d, in 
contrast to the scaling theory of localization, which per- 
mits only an insulator or a perfect metal ground state 
0. We then study the magnetoresistance and tunneling 
density of states in the presence of a magnetic field, and 
point out that these are excellent diagnostic tools to ex- 



tract key parameters and to test the applicability of the 
theory. At the end we shall discuss the MIT within the 
context of our theory of the metallic phase and comment 
on the effects of various symmetry breaking perturba- 
tions on the scenario we are proposing. Our main goal is 
to stimulate experimentalists to further study the metal- 
lic state both in the systems which have been studied up 
to now and possibly in other promising materials which 
we will discuss. 

We begin by summarizing the results of the scaling the- 
ory of interacting disordered systems jl0|-|l5| . In addition 
to the dimensionless resistance g, the theory is character- 
ized by the coupling constants 72, "f c and Z which obey 
the following scaling equations: 
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where y = — In A describes a rescaling of the length scale 
so that momenta in the range Afcp < k 2 < k 2 , are inte- 
grated out, where fco ~ (vpr) -1 is the short distance 
cut-off with t being the elastic scattering time. The 
parameter Z describes a rescaling of the energy scale, 
Z72 is related to the scattering amplitude in the triplet 
particle-hole channel, while Zj c is related to the singlet 
particle-particle (Cooper channel) amplitude. These pa- 
rameters can be interpreted in the context of Fermi liq- 
uid theory |l6|,[l7]]. For example, the specific heat linear 
T coefficient is modified by Z, so that Z plays the role 
of to* /m. The uniform magnetic susceptibility is given 
by Xs/Xs — ^(1+72) so that 72 plays the role of the 
Landau parameter — A° a . The key quantity in this the- 
ory is the diffusion propagation, which has a pole of the 
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form (Dq 2 — iZuj) -1 where D is related to the conduc- 
tivity a (which equals R^ 1 in 2d) by a = vqD, uq is 
the bare density of states. In the context of Fermi liq- 
uid theory, the diffusion pole can be written in the form 
(DQq 2 —ioj)^ 1 where Dq = D/Z has the interpretation of 
the quasiparticle diffusion constant. Equations (1-3) are 
derived to linear order in g and in the Cooper amplitude 
7 C but include all orders in the interaction amplitude 72. 
The exception is Eq. (4) for 7 C where the last term is 
quadratic in j c and independent of g. This term renor- 
malizes j c downwards, so that for j c > 0, 7 C becomes less 
important with scaling and can be neglegted for much of 
our subsequent discussions. The term 1 + 1 in Eq. (1) 
is written in a way to remind us that weak localization 
and singlet particle-hole channel in the case of Coulomb 
interaction give equal contributions to the enhancement 
of resistivity upon scaling. The next term is the contri- 
bution from the triplet particle-hole amplitude which has 
the opposite effect of reducing resistivity. According to 
Eqs. (2,3) both 72 and Z grow upon scaling. In fact, 
the growth is so rapid that they diverge at a finite scale 
yo, so that near yo they behave as 72 ~ (yo — y) -1 and 
Z — {Ho ~ y) 3 - This divergence signals the breakdown 
of the perturbative scaling equations. Here we want to 
make two important points: (1) the divergence of Z is in 
fact a necessary condition for the existence of a metallic 
state in 2d; and (2) due to the rapid growth of Z there 
is a wide range of temperature where the scaling equa- 
tions are valid and the system behaves like a metal. The 
key point is that the growth of Z forces us to perform 
scaling in an anisotropic manner in k space and energy 
space, a familiar situation in dynamical scaling. As we 
mentioned earlier, the key quantity is the diffusion pole 
(Dq 2 — iZui). The scaling procedure then consists of in- 
tegrating out the following regions in momentum space 
and energy space [^3| , 
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Xk 2 < k 2 < k 2 ; Afcp < —oj < k% . 

For Z growing with scaling, the energy or temperature 
scale decreases rapidly with scaling, and is given by 



T = XDk 2 /Z(X) 



(5) 



Strictly speaking, this formula needs further correction 
when Zi = Z(l + 72) becomes much greater than Z, 
because the energy denominator (Dq 2 — iZiuS) also ap- 
pears in some intermediate steps. However, the quali- 
tative point that the temperature scale can go all the 
way to zero remains valid. This is important because 
in one parameter scaling, the point has been made that 
the theory scales to either an insulator or a perfect metal 
(R n — > 0) in 2d, because the (3 function is always nonzero. 
The diverging Z at y = yo allows us to escape from 
this conclusion because in principle one can reach the 
point y = yo with g finite, so that according to Eq. (5) 
the system maintains a finite Rn as T — > 0. 



The next question is whether a metallic state can be 
realized in a region of parametric space and temperature 
where Eq. (1-3) are valid. From Eq. (2) and (3), it is 
apparent that the effective expansion parameter in the 
theory is (772. Then by starting with a sufficiently small 
g, it is possible to integrate Eqs. (1-3) until 372 becomes 
of order unity. Since Z diverges as (yo^y) 3 , much faster 
than 72 ~ (yo — y)~ l , the scaling can proceed to a rather 
low temperature before 1772 ~ 1 and the perturbative 
equations break down. By making the assumptions that 
g approaches a constant linearly in (y — yo) we conclude, 
using Eq. (5), that the low temperature behavior of the 
resistivity is given by R a (T) = R + cT 1 / 3 with c > 0. 
[Notice that at very low temperature, when 1772 ~ 1, 
the assumption that j c is negligible is no longer valid 
and indeed j c approaches a fixed point value 7* = 1 for 
72 — » 00. This would change the behavior of Z , leading 
to Z ~ (yo— y)~ J • This in turn modifies the temperature 
dependence of Rn = Ro + c'Ts when the regime 7 C ~ 1 is 
reached before getting out of the range of validity of Eqs. 
(1—4)]. To summarize, for sufficiently small g, we expect 
that initially g will exhibit InT correction over a broad 
temperature range. If 72 is sufficiently large to begin 
with, the InT correction is metallic-like. If 72 starts out 
small, the In correction resembles weak localization, but 
will change sign below a certain temperature scale when 
72 has grown sufficiently to overwhelm the localization 
term and the singlet contribution in Eq. (1). At a still 
lower temperature the resistivity drops rapidly, perhaps 
as T 1 / 3 (and possibly crossing over to Ts ) before the one 
loop scaling equation breaks down . This qualitative 
behavior has been confirmed [ p0[ by numerical integra- 
tion of Eq. (1-4). The point we wish to emphasize is that 
these equations predict a metallic behavior down to very 
low temperature in a region of parameter space where 
the one loop scaling equations remain reliable. Thus the 
existence of a metallic state over an experimentally ac- 
cessible temperture range should not in itself be a great 
surprise. 

We have seen that the key ingredient in arriving at a 
metallic state is the existence of a large 72. The question 
is whether 72 can be directly measured experimentally. 
We have mentioned that the uniform magnetic suscep- 
tibility provides a measurement of Z(l +72). However, 
this is a difficult, though not impossible experiment in 
a 2 dimensional electron gas (2lj. Instead, we find that 
magnetoresistance and tunneling in the presence of a par- 
allel field provide direct measurements of 72. A parallel 
field provides a Zeeman splitting of the spin states which 
cut off the S z = ±1 parts of the triplet particle- hole 
channel as well as the S z — part of the triplet and sin- 
glet particle-particle channel. This gives rise to positive 
magnetoresistance. The contribution coming from the 
particle-hole channel was calculated in the weak coupling 
limit in ref. |22|. This calculation was later extended to 
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strong scattering amplitudes 23 1 . Here we further extend 
this calculation to include the effect of the energy renor- 
malization Z . In analogy with Fermi liquid theory, we ex- 
pect the spin splitting of the quasiparticle to be given by 
£l s = (1 + 72)^ where ft s = ql^bH- Therefore the dif- 
fusion pole should be modified to (-Dq<7 2 — iuj — i&gSz) -1 
for the S z = ±1 components of the triplet particle hole 
channel. Inserting this modification into the expression 
for the S z = ±1 contribution to the conductivity, we find 



2 /.oo 



5a(T, H) = 



E 



, d ( w 
auj —r- I u> coth 



duo 



d 2 k 



D 2 n k 2 



2kTJ J (2tt) 2 q 
272 



~^ ±1 {D Q k 2 - ico- i£l s S z ) 2 D Q k 2 - i(l + j 2 )u - i^ s S z 

The parameters D, Z and 72 in this equation are scale 
dependent. Noting that the contributions for small H 
are dominated by small k and uj, we evaluate these pa- 
rameters at the scale A given by Eq. (5). The integrals 
are then performed following ref. p2f . In particular, we 
find that for small H, 
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We recover the weak coupling limit by setting 72 — > 
F/2 where F <C 1 is the interaction parameter in ref. 
p2[ . If we include the Cooper channel contribution, wc 
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will find an additional contribution of — 0.084^7 C (72 + 
l) 2 (gL/j,BH/kT) 2 . The above treats the effect of spin 
splitting only and is appropriate for H parallel to the 
plane. For perpendicular field we have, in addition to 
Eq. (7), the usual weak localization negative magnetore- 
sistance. In this case there is an additional contribution 
proportional to j c but now the orbital field scale given 
by fl h = 4:DeH/c also enters as a cut-off and the mag- 
netic field dependence from this term is more compli- 
cated. Since in the weak coupling regime we expect j c 
to scale to weak coupling, we shall concentrate on Eq. 
(7). The main point is that the quadratic in H term 
in parallel field magnetoresistance provides a measure- 
ment of the parameter 72. It will be very interesting 
to see if this parameter is indeed large in the metallic 
MOSFET samples and whether it increases with decreas- 
ing temperature. The available data are not systematic 
enough to answer these questions in the metallic regime. 
Most of the experiments on magnetoresistance are close 
to the MIT and for fields with Q s > kT. Qualitatively, 
the (positive) magnetoresistance increases as one moves 
away from the MIT @. This is in agreement with our 
expectation that 72 should consistently increase in order 
to establish a metallic phase. 

Another way to measure 72 is by tunnelling experi- 
ment. It was pointed out that the tunnelling density 



of states exhibit additional structure between the energy 
scales of the bare spin splitting ql^bH and the enhanced 
spin splitting due to interaction effects [p3[ . Following 
the Fermi liquid analogy, this second energy scale should 
be given by Q s . In particular, in 2d the derivative of the 
tunnelling density of states has logarithmic singularities 
at uj = gLfJ-BH and uj = {1 + J2)9lPbH . Thus tunnelling 
gives a direct measurement of 72. Recently a new tech- 
nique has been developed to tunnel into a 2d electron 
gas [p3 . It will be very interesting to apply it to the new 
metallic samples. 

As the field is increased, we expect a cross-over to the 
strong Zeeman splitting universality class. The detailed 
^cross-over is complicated, but the high field limit is one 
of the few fixed points which is controlled. The system 
always scales to an insulator, and in the weak disorder 
limit, a universal logarithmic temperature dependence 
was predicted @: a(T) = <7 + (e 2 /7r/i)(2-21n2) ln(Tr). 
As far as we know, this prediction has never been tested. 
The new MOSFET samples offer an ideal testing ground 
for this prediction. 

Up to now we have limited our discussion to the weak 
disorder case, when Eq. (1-4) remains valid. We now 
comment on the possibility of the existence of a non- 
trivial fixed point if somehow the scaling equations can 
be extended to strong coupling. In ref. [^9| the 2 loop 
contribution to the scaling equations was evaluated un- 
der the assumption of 72 3> 1 but for small gj2- The 
two loop scaling equations or ref. |l9| indeed exhibit a 
non-trivial fixed point. From this fixed point two separa- 
trices originate ending at 72 = and 72 = 00. Since the 
interesting part of the flow diagram is not in the weak 
coupling regime, the scaling equations and the details of 
the flow cannot be trusted. Nevertheless, the structure 
of the flow may be generic. Here we wish to make some 
general comments. If the initial 72 is not too large, the 
system exhibits a metal-to-insulator transition. An in- 
teresting feature of this flow is that on the metallic side 
of the separatrix the system reaches infinite 72 and Z at 
a finite scale A as in one loop order. Thus the discussion 
we gave earlier in this paper still holds and a metallic 
state with finite Rn is possible at T = 0. In fact, the 
metallic state in the low T limit exhibits a maximum 
metallic resistivity given by pm — {Trh/e 2 )gM, where gM 
is the value of g on the separatrix at 72 = 00. This 
gM is in general smaller than the value g* at the fixed 
point. Experimentally p* = (irh/e 2 )g* is determined as 
the resistance which separates the metallic and insulat- 
ing states at higher temperature. This feature seems to 
be consistent with currently available data. For example, 
the data of ref. § yields p M » 0.1^ and p* w 2^. 

The scaling behavior near the MIT will be controlled 
both by the existence of a fixed point at finite g* and 
72 and by the runaway towards g ~ gM and 72 = 00. 
Then one can show that Rn = p(T/(Sn)" z ) where 5n is 
the deviation from the critical density and the critical 
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indices v and z are determined by the fixed point, p is a 
scaling function and according to the previous discussion 
p(oo) = {ith/e 2 )g* and p(0) = (irh/e 2 )g M - 

Besides the magnetic field, other symmetry breaking 
perturbations have relevant effects on our picture of the 
2d metallic phase. Spin flip scattering by magnetic impu- 
rities will cause a crossover to a low T insulating phase. 
The effect of spin orbit (SO) scattering is more intrigu- 
ing. In d — 2, intrinsic SO coupling or SO scattering by 
impurities only affects the out of plane component of the 
spin p5| . In this case the one loop equations ]26| still 
lead to a diverging behavior of the (S z — 0) triplet am- 
plitude and a metallic phase at low T. We suggest that 
the above discussion on the MIT applies in this case even 
though the 2d SO could result into a different universal- 
ity class. A much more dramatic effect on our theory 
of the metallic phase is the SO scattering deriving from 
possible asymmetry of the confining potential since it is 
equivalent to a 3d SO coupling and cutoff all triplets [^7) . 
If this coupling is sizeable, the theory predicts an insu- 
lating behavior at zero temperature |p8| , at least in the 
limit in which the SO band splitting is less than the in- 
verse elastic scattering time. In our opinion, evidences of 
2d or 3d SO are still lacking. 

The scenario wc outlined in this paper has the advan- 
tage of permitting a metallic state in 2d and therefore a 
metal-insulator transition. However, given the uncertain- 
ties of the strong coupling theory, a good strategy is to 
approach the MIT from the metallic side and try to gain 
a thorough understanding of the metallic state. This mo- 
tivates us to propose magnetic susceptibility, magnetore- 
sistance and tunnelling experiments as ways to directly 
measure the key parameters of the theory 72 and Z . We 
also worked out the qualitative behavior of the tempera- 
ture dependence of the resistivity, in a regime where the 
theory is valid. Here our results do not compare favorably 
with experiments. The data of ref. [0 and ref. jl| have 
been fitted to the form p(T) = p + exp(-T /T). This 
is very different from the InT dependence followed by 
a low temperature power law that we predict. Further- 
more, the parameter T appears to scale with the Fermi 
energy which is relatively small in these low density sys- 
tems. Thus the possibility remains that some physics on 
the scale of the Fermi energy is playing the dominant role 
and the data are far from the low energy scaling regime 
we considered here. We believe these questions can be 
addressed by more detailed studies of the metallic state 
along the lines suggested in this paper. Yet another pos- 
sible research direction to confirm the theory here pre- 
sented is to study 2d systems where 72 is expected to be 
large to begin with, such as almost ferromagnetic metal- 
lic thin films. Examples are weak ferromagnets such as 
MnSi or TiBe2, if the ferromagnetism can be suppressed 
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